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We derive the exact generating function for planar maps (genus zero fatgraphs) with ver- 
tices of arbitrary even valence and with two marked points at a fixed geodesic distance. 
This is done in a purely combinatorial way based on a bijection with decorated trees, lead- 
ing to a recursion relation on the geodesic distance. The latter is solved exactly in terms of 
discrete soliton-like expressions, suggesting an underlying integrable structure. We extract 
from this solution the fractal dimensions at the various (multi)-critical points, as well as 
the precise scaling forms of the continuum two-point functions and the probability distribu- 
tions for the geodesic distance in (multi)-critical random surfaces. The two-point functions 
are shown to obey differential equations involving the residues of the KdV hierarchy. 
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1. Introduction 

The study of the statistical properties of random graphs is relevant for many problems 
in physics, such as two-dimensional quantum gravity or fluid membrane statistics. By 
random graphs, we here mean graphs embedded in a surface of given genus, also known 
as fatgraphs or maps. Indeed, these are the natural discretization of fluctuating surfaces, 
on which matter systems may be defined. Much is known to this day on the enumeration 
of maps of fixed genus either by combinatorial techniques [|I] or by the use of matrix 
integrals [0-^]. When applied to graphs of large size, these results give rise to various 
scaling behaviors depending on the critical universality class at hand, described in the 
continuum by the coupling to 2D quantum gravity of some 2D Conformal Field Theories 
(CFT) with central charges c < 1. In particular, the critical behaviors are characterized by 
the famous KPZ scaling relations ||. As an example, the so-called one-matrix model, which 
enumerates fatgraphs with vertices of arbitrary valence and of any fixed topology, displays 
a set of multicritical points corresponding to the CFT's with central charges c(2, 2m + 1), 
m = 1,2,3,.... The first value m = 1 corresponds to c = 0, i.e. a model of pure 2D 
quantum gravity without matter, representing the universality class of generic random 
maps of large size. The higher order multicritical points may be reached by weighting 
vertices according to their valence, and by suitably fine-tuning the vertex weights. 

Most results obtained so far concern global properties of the graphs which, in the 
continuum, can be translated into correlation functions integrated over the positions of 
their insertion points on the surfaces. Little is known however on more refined properties 
of random graphs such as the dependence of correlators on the distances between their 
insertion points. By distance we mean here geodesic distance along the graph, namely the 
length of any shortest path between say two given faces. These refined properties seem to 
remain beyond the reach of the standard matrix model treatment. Still, some results were 
obtained in Refs.|6|-|8|], where in particular combinatorial arguments were used to derive 
the universal scaling two-point function for critical genus zero surfaces with two marked 
points at a fixed geodesic distance. However, an explicit form for the two-point function 
was obtained only in the pure gravity case and close to the critical point. 

More recently, using a completely different combinatorial approach, the same result 
was recovered in Ref. |§ in the case of random tetravalent planar maps. This work relies on 
the existence of a bijection between planar tetravalent maps and rooted trees with vertices 
labelled by the geodesic distance to the root. The construction however is specific to the 
tetravalent case, and does not provide a closed form for the discrete solution either. 
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In this note, we address the general question of the enumeration of planar maps with 
two marked faces at a fixed geodesic distance. In the case of maps with inner vertices of 
arbitrary even valence, we derive explicit expressions for the generating function G n of 
"two-leg diagrams", namely planar maps with two distinguished univalent vertices (legs) 
whose adjacent faces lie at a geodesic distance n of one-another, and with weights gi per 
inner 2z-valent vertex, % = 1,2,3,.... These results are obtained via the use of a more 



general bijection between two-leg diagrams and decorated trees first found in Ref. [10 1 and 



extended in Ref. [|TT[ , which we adapt so as to keep track of the geodesic distance between 
the two legs. This leads to a simple algebraic recursion relation on n, which allows to derive 
a closed expression for the solution G n . We give several explicit expressions depending 
on the maximum valence of the vertices of the graph. With these solutions in hand, 
we easily derive the continuum scaling two-point functions corresponding to the various 
multicritical points. This allows to recover the results of Refs.[[7] and || in the generic 
critical case, and also provides the generalization to all higher order multicritical points. 
We also extract critical exponents such as the fractal dimension in various approaches 
to the multicritical points, as well as probability distributions for the (rescaled) geodesic 
distance in multicritical planar maps of fixed but large size. 

The paper is organized as follows. In Sect. 2 we detail the general correspondence 
between two-leg diagrams and decorated trees. Sect. 3 is devoted to the incorporation of 
the geodesic distance n between the two legs in this setting and leads to a general recursion 
relation on n. This recursion relation is solved first in the tetravalent case (Sect.4.1), and 
exploited to derive a critical fractal dimension (If = 4 (Sect. 4. 2) and to recover the known 
continuum two-point function for critical surfaces (Sect. 4. 3). We also use our explicit 
solution to derive in Sect. 4. 4 the probability distribution for rescaled geodesic distances 
in two-leg diagrams of fixed but large size. The general case of maps with valences up to 
2(m + 1) is solved in Sect. 5, where we discover a remarkable connection with m-soliton 
solutions of the KP hierarchy. This general solution is applied in Sect. 6 to the so-called 
"hard dimer" model, which displays the first higher order critical point (tricritical point). 
In Sect. 6.1 we detail the solution and display in particular the first few G„'s. This leads 
to the computation in Sect. 6. 2 of the fractal dimension dp = 6 and the continuum two- 
point function for tricritical surfaces, corresponding to a generic approach to the tricritical 
point. We also discuss the case of a non-generic approach along the critical line, which 
displays another fractal dimension dp = 4 and a different two-point scaling function. 
We finally obtain in Sect. 6. 3 the probability distribution for rescaled geodesic distances 
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in tricritical surfaces of fixed but large size. All these results are extended to the case 
of the m-th order multicritical point in Sect. 7, where we obtain the fractal dimension 
cIf = 2(m+l) (Sect. 7.1), the generic continuum multicritical two-point function (Sect. 7. 2) 
and the corresponding probability distribution (Sect. 7. 3). We gather a few concluding 
remarks in Sect. 8, namely the relation to matrix models (Sect. 8.1), the extension to two- 
point functions of other (less relevant) operators (Sect. 8. 2), some generalizations to other 
classes of planar graphs such as constellations (Sect. 8. 3), as well as a connection to the 
so-called Integrated SuperBrownian Excursion (ISE) (Sect. 8. 4). We also gather a few 
technical details in appendices A and B. 



2. Planar maps and decorated trees 

There is a general correspondence between planar maps and special classes of deco- 



rated trees ]TQ,|TTj] . In this paper, we concentrate on the simpler case of maps with only 
vertices of even valence. In this case, there is a bijection between so-called "two-leg di- 
agrams" and so-called rooted "blossom" trees as defined below. By two-leg diagram, we 
mean a planar map with two external legs (i.e. two extra univalent vertices) distinguished 
as in- and out-coming, whereas all inner vertices have even valences. As opposed to the 



definition used in Ref. |Tl|, we do not require that the two external legs lie in the same 
face. By convention, we decide to always represent the diagram in such a way that the 
in-coming leg is adjacent to the external face in the plane. On the other hand, by rooted 
blossom trees we mean rooted planar trees with two kinds of endpoints called buds and 
leaves, with only inner vertices of even valences and such that each inner 2/c-valent vertex 
is connected to exactly k — 1 buds. It is easily seen that any such tree has exactly one 
more leaf than bud. 

The one-to-one correspondence is illustrated in Fig.l. To each two-leg diagram we 
associate a rooted blossom tree as follows. Starting from the in-coming leg, we iteratively 
visit all edges of the diagram adjacent to the external face in counterclockwise direction. 
At each step we cut the edge iff the remaining diagram remains connected. The first half 
of the edge is replaced by a bud, the second one by a leaf. This procedure is repeated 
with the remaining diagram until a tree is obtained. Finally, the in-coming leg is replaced 
by a leaf, while the out-coming one is chosen to be the root. It is shown in Refs. ]TI1| , P^T| 
that the resulting tree is indeed a blossom tree as defined above. This transformation 
is moreover invertible, by iteratively connecting each bud to the closest available leaf in 



3 



(a) 



(b) 



III/ 




III/ 




(c) 



(d) 



Fig.l: Illustration of the one-to-one correspondence between two-leg dia- 
grams and blossom trees. Starting from a two-leg diagram (a), we apply the 
iterative cutting procedure explained in the text, which here requires two 
turns around the graph. The edges cut during the 1st and 2nd turn are dis- 
played in (b) with the corresponding index 1,2. Each cut edge is replaced 
by a bud-leaf pair (c), while the in-coming leg is replaced by a leaf and the 
out-coming one by a root, finally leading to a blossom tree (d). Conversely, 
the matching of buds and leaves of the blossom tree (d) rebuilds the edges of 
(a). 



counterclockwise direction and fusing these bud-leaf pairs into edges. This construction 
leaves one unmatched leaf, adjacent to the external face, and that we replace by the in- 
coming leg, while the root becomes the out-coming one. Note that in the above fusion 
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procedure the root may possibly be encompassed by a number of nested edges, in which 
case it does not lie in the external face. The number of "encompassing" edges will be called 
the depth of the root in the blossom tree. 

The generating function R for blossom trees with a weight gk per 2/c-valent vertex 
satisfies 



as obtained by inspecting all possible configurations of the vertex attached to the root. 
The function R is uniquely determined as the formal power series solution of eq.( [2.1|) in 
the gffc's with constant term 1. By the above correspondence, R also coincides with the 
generating function T^i of two-leg diagrams with a weight gk per inner 2/c-valent vertex. 
This is to be contrasted with the generating function 1?2 for two-leg diagrams with both 
legs in the external face, or equivalently blossom trees with root at depth 0, which reads 



3. Geodesic distance 

A nice feature of the above bijection is that it keeps track of the geodesic distance 
between the in- and out-coming legs in the two-leg diagrams, defined as the minimal 
number of edges crossed by a curve connecting the two legs (for instance, the two legs of 
the diagram of Fig.l (a) are at distance 1). Indeed, this distance is nothing but the depth 
of the root in the corresponding blossom trees. More precisely, it was shown in Refs. [|10| - |T2| 
that all the edges cut in the cutting algorithm belong to minimal paths (i.e. paths crossing 
a minimal number of edges) going from the external face to all internal ones. In particular, 
the path from the external face (adjacent to the incoming leg) to that adjacent to the 
out-coming leg and crossing each encompassing edge exactly once is minimal. 

This leads to the refined definition of the generating function G n for blossom trees 
with root at depth n, which alternatively generates the two-leg diagrams with geodesic 
distance n between the legs. In particular, we have 1?2 = Go and I\i = ^ n>0 G n . In the 
following it will prove easier to use the generating function R n for blossom trees with root 
at depth at most n, such that 




(2.1) 
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(2.2) 



G 



n 



Rn — R. 



■n— 1 



(3.1) 



IYi = R = lim R 



n 



n — >oo 
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with the convention that _R_i = 0. 

The first result of this paper is a refined version of eq.( ^.l| ) which keeps track of the 
depth of the root of the trees. 



ILL/ 





(a) 



(b) 




(c) 



Fig.2: The contour walk (b) associated to a blossom tree (a) obtained by 
visiting the buds and leaves around the tree in clockwise direction starting 
from the root. Each bud (resp. leaf) corresponds to a +1 (resp. —1) step. 
The walk is decomposed in (c) according to the subtrees attached to the root 
vertex, here two buds and three blossom trees of respective root depths 0, 
2, 0, viewed as the maximal relative heights reached by the corresponding 
portions of walk. Here, the depth of the root is 3, which corresponds to the 
maximal height attained by the contour walk. 



For any blossom tree, let us visit all its buds and leaves successively in clockwise order 
around the external face starting from the root. This defines the contour walk of the tree, 
namely a walk on the integer line starting from and moving up (height shifted by +1) 
for each bud encountered and down (height shifted by —1) for each leaf. As leaves are in 
excess of one on each of these trees, the contour walk stops at height —1. We now easily 
identify the depth of the root as the maximal height reached by the walk (see Fig.2). 

As before, a recursion relation is obtained by inspection of all possible environments 
of the root in a given blossom tree. In the case of the smallest blossom tree made of a 
single leaf, the depth of the root is 0, hence this tree contributes 1 to R n for all n > 0. 
In the case of a tree with a 2/c-valent vertex attached to the root, we have exactly k — 1 
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Ill/ 




Fig.3: The root vertex environment is characterized by the cyclic sequences 
of attached buds and blossom subtrees, represented with a circle. The index 
mi denotes the number of buds immediately before the i-th blossom subtree, 
whose root depth is denoted by pi. The corresponding countour walk is also 
sketched, and gives a visual representation of eq.(3.2). 

buds and k descendent blossom trees attached (see Fig.3). Turning around the vertex 
in clockwise direction from the root, let m,j denote the (possibly vanishing) number of 
consecutive buds immediately before the z-th descendent blossom tree, i = 1, 2, k (with 
necessarily ^2mi < k — 1, rrii > 0). Let moreover pi, i = 1,2, ...,k denote the depth of 
the root of the i-th descendent blossom tree, considered as an independent rooted blossom 
tree. With these definitions, the depth p of the root of the whole tree reads 



This relation is clear in the contour walk picture, where the walk associated to the whole 
tree decomposes into +1 steps for the buds attached to the vertex, and into subwalks 
associated with the descendent blossom trees, each contributing by a global — 1 height 
shift. The quantity pi + rrii + rrii-i + ■ — h mi — (i — 1) is nothing but the maximal height 
reached within the i-th subwalk, which starts from height m,+mi_i+- • -+mi — (i— 1). From 



p = max {pi + rrii + m,_i H h mi - (i - 1)} 



(3.2) 
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eq. ( |3.2|) , we see that the condition p < n is equivalent to Pi < n + i — l — mi— mi — • • • — 
for all i = 1, 2, A;, which finally translates into the recursion relation 



A- 



n+i-l-mi mi 

(3.3) 



k>l m,i>0 i=\ 

m^ + ...+m^</c — 1 



with the convention R n = for n < 0. This may be viewed as a generalization of eq.( |2.1| ) 
which is recovered in the limit n — > oo, where _R n — ► i?. 

The equation ( |3.3|) may be rewritten in an algebraic way by use of a shift operator a 
acting on the canonical basis {p n }nez for sequences as ap n = p n +i, and of its inverse a -1 
such that o~~ 1 p n = p n -i, namely introducing the operator 

Q = a + a~ 1 r (3.4) 

where r acts diagonally as rp n = R n p n , eq. (|3.3[ ) takes the form 

Q n - hn = l + J29k(Q 2k ~ 1 )n-l,n (3.5) 
k>l 

where Qij stand for the matrix elements in the canonical basis {p n }- This algebraic 



formulation (|3.5| ) is very similar to the main recursion relation obtained in the context of 



one-matrix models with M x jV Hermitian matrices and with even potentials of the form 
V(M) = H — J2k>i 9 k ^ik~i known to generate fatgraphs with 2/c-valent vertices weighted 
by Qk 0- In this setting, the operator Q represents the multiplication by an eigenvalue A, 
acting on polynomials of A. In the basis made of the monic polynomials p n (A) orthogonal 
with respect to the measure d\e- Mvw on IR, Q takes the form and the recursion 
relation for the i? n 's takes precisely the form ( |3.5| ), except that the first term on the rhs 
is changed from 1 to nj M . This coincidence is both remarkable and mysterious as in the 
matrix model recursion, the i? n 's give access to the generating functions for higher genus 
graphs, as opposed to the present ones, which only concern planar graphs, but on the other 
hand allow to explore geodesic distances, a task out of reach in the matrix model. 



4. Tetravalent case 

4-1. Exact solution 

Let us now show how eq. (jO)) may be used to extract exact expressions for the i? n 's. 
As a first simpler but instructive case, let us restrict our study in this section to the case 
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of tetravalent maps, i.e. two-leg diagrams having only tetravalent inner vertices. Going to 
generating functions, this amounts to taking g^ = gSk,2, in which case eq. (pTT|) reduces to: 



R = l + 3gR 2 (4.1) 

with solution 

6g 

which displays the well-known critical value g c = 1/12 of g for pure quadrangulations. 
More generally, the relation (|3.3|) reduces to: 



R n — 1 + gR n (R n -i + R n + R n +i) (4-3) 

Conceptually, this equation, valid for all n > 0, may be viewed as a recursion of order 2 in 
n, fully determining all the -R n 's from the initial conditions: 

(i) R-i = 0; 

(ii) Rq = ]?2, where 1?2 is obtained from eqs. (|2.2| ) and ( |2.1|) , namely 1?2 = R — gR s with 
R as in (O). 



From the combinatorial origin of the equation, we expect that the non trivial initial con- 
dition (ii) may be replaced by the alternative boundary condition at infinity that: 
(ii') lim n ^ 00 R n = R with R as in (|42|) . 



The general solution to eq. (|4.3|) satisfying (ii') reads: 



- p (l-Ax n+1 )(l-Ax n + 4 ) 
n " (l-Ax-+ 2 )(l-Ax-+3) ( > 

where A is an arbitrary "integration" constant and x is the solution of the characteristic 
equation: 

= x + ~ 4 - 5 

gR x 

with |x| < 1. Note that for < g < g c , the lhs is larger or equal to 2, hence x is real and 
positive. The condition < x < 1 (which ensures the convergence (ii') of R n to R) can 
always be satisfied for g < g c . 

That the expression ( |4.4[ ) solves eq. ( |4.3| ) may be checked explicitly. In a more construc- 
tive way, this solution is obtained by expanding R n in the vicinity of R as R n = R(l — p n ) 
and first expressing eq.( [4.3| ) at first order in p n with the resulting linearized equation 

p« = gR{p^_ x + ti ] + til,) + 3gRpV (4.6) 
9 



The general solution of this linear recursion such that lim n ^ 00 p£ = has the form 

Pn = ot\x n where x is the solution of the characteristic equation (|4.5| ) with |x| < 1. 
Further expanding R ni we have 

Pn = J2^(x n y (4.7) 

i>l 

Substituting R n = R(l — p n ) and picking the term of order k + 1 in x n yields the relation 



a k+1 l x k+1 



)k — 1 -. u 
1 — X 
- 4.9 
1 — x 



solved recursively as 



Picking a i = x(l — x)(l — x 2 )A, R n is easily resummed into (|4.4j 
Further imposing the initial condition (i) fixes A = 1, hence 

[1 - x n+1 )(l - x n+4 ) 



Rn — -R- 

= R- 



1 - x n+2 )(l - x n + 3 ) 
U n (w)U n+3 (w) 



(4.10) 



t/" n+ i(ty)t/" n+2 (w) 

where the U^s are the Chebyshev polynomials of the first kind, namely obeying U n +i(w) 
wU n (w) — U n -i(w) and Uq = 1, U\ = w, expressed in the variable 



A first outcome of this solution is that it allows to recover the value of T2, namely 

r 2 = Ro = Rt^tt = *^rf = R - gR s (4.12) 
1 + x + | 1 - 3gR 

as expected, which a posteriori shows the equivalence (ii) <-> (ii'). 

As an illustration of the result ( |4.10| ), we list below the first few generating functions 
G n (g) = R n — R n -i for two-leg diagrams with legs at a geodesic distance n = 0, 1, 2... 



G Q (g) = 24 G 9 n 1+ ^LJ ^- = l + 2g + 9g 2 + 54g 3 + 37% 4 + 2916<? 5 + • - ■ 
90(1 + V 1 - 12 #) 

27o-2 + (2-15o)v/l-12o , , . . 

Gi(g) = |= j, \ - J^L^, =g + Sg 2 + 65g 3 + 554/ + 4922^ + • • ■ 

27g(l -8g + V 1 - 12 #) (4.13) 

_ 252^ 2 - 6^ - 1 + (1 + 12^) V ^T2^ 

2 ^ J 27(190^ 2 - 510 + 3 + (3 - 33g + Ug 2 )y/1 - 12g) 

= g 2 + 15g 3 + 179# 4 + 1995# 5 + • • • 
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Fig.4: The tetravalent two-leg diagrams with two inner vertices, classified 
according to the geodesic distance n between the two legs, namely n = (a), 
n = 1 (b) and n = 2 (c). 

We have represented in Fig.4 all the two-leg diagrams with 2 vertices, among which 9 (resp. 
8, 1) have the two legs at geodesic distance (resp. 1, 2), in agreement with the g 2 terms 
in eq.(|7L3|). 

4-2. Fractal dimension 



As a direct application of the solution Q4.10| ), let us evaluate the "fractal dimension" 
of large tetravalent planar maps, defined via the number R Uj n of two-leg diagrams with 
distance less or equal to n between the two legs and with iV vertices. Indeed, at large N, 
the limiting ratio 

Bn = lim (4.14) 



lim ^?> N 

11 



may be taken as a good estimate of the average number of points at a geodesic distance 
less or equal to n from a given point in random tetravalent maps of infinite size. It is 
expected to behave like 

B n ~ n dp for large n (4-15) 

where dp is the fractal dimension. 
We have 

Rn > N = j 2^F+T^ (416) 
where we must use the expression ( PH ) for R n , in which we substitute expression ( |4.2|) for 
R(g) and solve eq.( ^.5| ) for x as a function of g. Upon performing the change of variable 
g — > v = gR(g), namely g(v) = v{l — 3v), we arrive at 

r dv(i-6v) i (i- X (vr+ i )(i-x(vr+*) 

n ' N J 2in(v(l-3v)) N + 1 l-3v(l-x(v) n + 2 )(l-x(v) n + 3 ) { ' 

where we have used R(g(v)) = 1/(1 — 3v) and the expression x = x(y) = (1 — 4v — 
V^-^v~^-12v^) I (2f ). The large iV behavior is obtained by a saddle-point approximation, 
with the result 

(12) N (n + l)(n + 4) 

Rn n ~ const. {) {(140 + 270n + 179n 2 + 50n 3 + 5n 4 ) (4.18) 

iV"! (n + 2)(n + 3) y 

which finally gives the ratio 

280 (n + 2)(n + 6) 56 
which yields dp = 4, as expected for pure gravity. 

4-3. Critical scaling 

A continuum limit may be reached by letting g tend to its critical value g c = 1/12, 
corresponding to gR — >• 1/6 and x — > 1. More precisely, we write 

<7 = ^(l-e 4 ) 9 R=\{l-e 2 ) (4.20) 

from eq.(^.2j). In turn, the characteristic equation ( fO| ) yields 

x = e - ae + 0(e 2 ) a = V6 (4.21) 
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As seen from eq.( f4.10 ), a sensible limit is obtained by writing 



r 



n=- (4.22) 



e 



and letting e — > 0. Writing the scaling variable r as r = n/£, we see that e plays the 
role of the inverse of the correlation length £. As we approach the critical point, we have 
£ = e _1 = ((<7c — g)/g c ) V with a critical exponent f = 1/4, in agreement with v = 1/d.F, 
as expected from general principles. Performing this limit explicitly on the solution ( |4.10| ) 
yields an explicit formula for the continuum partition function T(r) of surfaces with two 
marked points at a geodesic distance larger or equal to r: 



>6 e 2 R ~~ ~ d^ 2 yw ' z ' ' sinh 2 (^f r) 



T{r) = lim — = -2-^Log sinh (^fr) = , — (4.23) 



Upon differentiating with respect to r, we obtain the continuum partition function for 
surfaces with two marked points at a geodesic distance equal to r: 

g(r) = -f{r) = 3V6 ^fl^ (4.24) 

sinh (vfrj 

in agreement with the form obtained in Ref. @ . 

Note that the precise form of the scaling function J-{r) may alternatively be obtained 
by solving the continuum counterpart of eq.( f4.3|) . Indeed, writing 



R n = R(l -t 2 T{ne)) (4.25) 
and expanding eq. (|4.3|) up to order 4 in e, we obtain the following differential equation 

F"{r) - 3^ 2 (r) - 6^(r) = (4.26) 



It is easy to check that F(r) as given by ( fl.23|) is the unique solution of (|4.26| ) with 
boundary conditions J-(r) — ■> oo when r — > and J-{r) — > when r — > oo. 

4-4- Continuum probability distribution for geodesic distances 

The scaling function Q (r) obtained in the previous subsection is the continuum "two- 
point function" for surfaces with two marked points at geodesic distance r. It is obtained 
in a grand-canonical formalism, in which the area of the surfaces (number of vertices 
in the corresponding maps) is not fixed. In particular, Q(r) cannot be interpreted as 
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a probability distribution as it is not normalizable. Such a problem does not occur in 
the canonical formalism, i.e. when working with maps of fixed area N. A well defined 
probabilty is given by R n ,N / Roo,n, where Rn,N is the number of two-leg diagrams with A 
inner vertices and legs distant by at most n, as before, and Roo,n is the total number of 
two-leg diagrams with iV inner vertices. Starting from the expression Q4.16Q for R n) N and 
performing again the change of variables g = v(l — 3i>), we get 



*-» = # K^-S)™ **™ (4 ' 27) 



where the integral is to be taken on a small contour around the origin. At large N, the 
integral is dominated by the saddle-point v = V c = 1/6. We therefore perform the change 
of variables v = V c (l + ut/vA), with u to be integrated on the real line when A — ► oo. As 
opposed to the computation of Sect. 4. 2 where R n ,N was estimated for finite n and large A, 
we may obtain a non-trivial distribution by letting A — > oo while n scales as n = aA 1 / 4 , 
in which case we may use eq. fl4.25p with e = yf—iu/N 1 ^, namely 



R n (v) =r(i + i-j=F(ay/^iu) + 0(i)J 

= I i 1 + W + °&) + 'TN^aV^-u) + 0(1,) (4.28, 

= 2 ( 1 + ^( 1 + ^ (Q ^ , )) +0( ^ ) 
Substituting this into eq. ( |4.27|) , we find 

Rn,N ~ 2 ^7^ J duu 2 e- u (l + F{aV^iu)) 



- oo 



7TAT3/5 



= \— — : J duu 2 e- u (l + Re \F(a\/^iu. 



(4.29) 



and a similar expression for the limit Roo,Ni obtained by taking T — > in the above. 
This leads to the probability $(a) for a random two- leg diagram to have legs at a rescaled 
geodesic distance less than a: 

"°° 2 _ u 2 cosh(2ay / 3w) + cos(2ay / 3w) + 8cosh(av^3w) cos(a\/3u) — 10 

Sir \ Ql J — J / CLUjU G 2 

(cosh(av / 3w) — cos(av / 3w)) 

(4.30) 
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Some remarks are in order concerning this explicit form. First we may explicitly expand 
the probability distribution $(a) as a formal power series of a, with the result 

*(a) = sf^(-l) j+1 a Aj S 2j -^B 4j+2 (4.31) 

where Bj denote the Bernoulli numbers. This series has a radius of convergence zero 
but provides an asymptotic expansion valid to any finite order. In particular, the lead- 
ing small a behavior is ~ 3a 4 /28, which translates into the equivalent R n ,N ~ 
(3/28)(n 4: /N)R 00j N for finite n and large N, while from eq.( [4.19| ) we have R n ,N ~ 
(3/56)n 4 i?o,Af • These relations are compatible, as we know that i?o,iv = 2i2 00) jv/(iV + 2) 
from combinatorial arguments ||1U|| . 

Secondly, upon differentiating with respect to a, we obtain the probability distribution 
p(a) for two-leg diagrams with legs at rescaled geodesic distance a, with the result: 

4 f°° 2 

p(a) = / duu 2 {2u 2 - 3)e~ u2 x 

V 71 "" Jo 

cosh(2av / 3w) + cos(2av / 3w) + 8 cosh(av / 3w) cos(av / 3w) — 10 

X ', " — 2 

(cosh(ay / 3w) — cos(a\/3u)) 

(4.32) 

At large a, the function p(a) decays as exp(— Ca s ), with C = 3(3/8) 2 / 3 and with a critical 
exponent 6 = 4/3, in agreement with Fisher's law S = 1/(1 — v), a result already observed 
in Ref.0. 

O(oc) p(oc) 




(a) (b) 

Fig.5: Plots of the probability distributions $(ct) (a) and p(a) (b) for two-leg 
diagrams with legs at rescaled geodesic distance at most (a) or equal to (b) 
a. 
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The distributions and p(a) are plotted in Fig. 5 for illustration. 



5. Exact solution in the general case 



Let us now turn to the general case of maps with bounded valences, namely with 
arbitrary even valences 2, 4, 6, ... up to say 2(m+ 1). This corresponds to keeping arbitrary 
vertex weights gi,g%, ...,g m +i while setting gj = for j > m + 1. The equation ( |3~"3l ) 
becomes a recursion relation of order 2m namely expressing R n +m in terms of R n +j, 
j = m — 1, m — 2, — m. This equation is valid for all n > 0, upon taking the initial 
condition 

(i) = R-2 = ••• = R-m = 0. 

To completely fix the solution, the condition (i) must be supplemented by 

(ii) the data of the m first terms Ro, R±, R m -i- 

The explicit combinatorial expression for Rq is known from the identification with the 
two-leg function 1?2, given by eq.( |2.2[ ). As to the other Rj, j = 1,2, ...,m — 1, no direct 
combinatorial expression is available yet. In view of the previous section, and from the 
combinatorial interpretation of R n , we replace the initial data (ii) above with the condition 
(ii') lim n ^ oc R n = R, where R is the solution of eq.(|2~J]). 

As we shall see, this latter condition is indeed sufficient, together with initial condition (i) 
to completely fix the solution R n . 

To reach this solution, we follow the same strategy as in the tetravalent case of Sect. 4. 
We first obtain constructively the solutions to eq. (|3.3|) with the only boundary condition 
(ii'), expressed in terms of m free "integration constants" Ai, A m . We finally implement 
the boundary condition (i) to further fix these constants. 

Starting with the condition (ii'), let us first expand R n around its limit as R n = 



R(l — p n ) and express eq. (|3.3| ) at first order in p 7 



P (1) 

rn 



E. 



,9k+iR k 2J 

fc=0 j=-k 




(5.1) 



k=0 j=-k 

The coefficient of p n X +j lIV the first line is easily obtained in the contour walk picture as 
the number of walks from height n to height n — 1 with 2k + 1 steps and a marked step 
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down from height n + j, itself decomposed into a walk from height n to n + j and say I 
steps, and a walk from height n + j — 1 to n — 1 with 2k — I steps. Looking for a solution 
Pn+j = a i xU i we get the characteristic equation 

J^f2k + 1\ 1 1 - x 2k+1 ~ 21 
Xm(x)=0, Xm(x) = 1 -22g k+1 R k 22 ( , J^fcZT 

I — n ; — n \ / 

(5.2) 



fc=o ;=o 

m k 



i — n 1— n V / 



fc=0 Z=0 

where the IPs are the Chebyshev polynomials of the first kind, with again w = ^/x + l/^/x. 
Note that Xm( x ) is invariant under the interchange x «-> and moreover is a polynomial 
of degree m in x + 1/x. The most general solution p^ is a linear combination 

m 

p<P =Y,<*i*i (5-3) 

where j = 1, 2, m are the m zeros of % m with modulus less than 1. Further expanding 
p n order by order in the x], as p n = E Ill ,...,n ra >oS,.,« m W) ni -(C)''"' wi th «o,o,...,o = 
and ao,..,o,i,o,...,o = a j f° r the 1 in position j, we get a recursion relation for the coefficients 
«ni,...,n m , which determines them all in terms of the initial integration constants ai, a m 
of ( |5.3|) . Upon redefining them like in previous section as ai = X{(1 — Xi)(l — x?)\i, and 



resumming the expression for R n , we arrive at 



(m) (m) 
p p Un U n+3 

n ' n ^ (m) (m) 



U n+l U n+2 



i m) =E(- i ) z e ri A -< m n (5 - 4) 

;=0 Kmi<...<mi<m i=l l<*<i<' 



(1 X a X}^j 



where, as before, the must be chosen so that Xm( x i) = 0, \xi\ < 1, % = 1,2, ...,m. 
Proving this result is beyond the scope of the present paper. That the form ( |5.4|) actually 
solves ( |3.3|) with x's solutions of eq. (|5.2|) with modulus less than 1 may be explicitly checked 
for small values of m = 2, 3, 4. 

A few remarks are in order regarding the particular form (|5.4| ). First the expression 
as a function of R and of the Xj's is completely independent of the particular values of 
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the parameters g^. The link with the specific eq.( [3.3|) appears only through the expression 
for R solving ( |2.1] ) and the characteristic equation (|5.2| ) satisfied by the x's. The form of 



is reminiscent of that of the so-called N-soliton tau-function of the KP hierarchy [JT3| 
with the generic form 



C i,j 



(5.5) 




(Pt ~Qj){Qi ~Pj) 



where the 77's contain the times' dependence of the KP hierarchy. We see that, at least from 
a formal point of view, our solution (5^4) corresponds to taking N = m, e Vi = — x™ +m A 



l ! 



Pi = Xi and qi = 1/xt. This relation, though suggestive, remains mysterious at present. 
Secondly, we note that any "N-soliton" solution of the form ( |5.4| ) with m — > N > m reduces 
to the m-soliton solution when we pick the x's among the m zeros of Xm with modulus 
less than 1. Indeed, for N > m, some of the x's will be equal resulting in the vanishing 
of the corresponding c^j's. Analogously, taking an JV < m solution is perfectly licit, as 
it amounts to taking some A's to zero, but it corresponds to a less generic solution. In 
particular, the one-soliton solution ( |4.4|) is solution to all general equations, provided x 
solves the corresponding characteristic equation. Thirdly, inspired by the determinantal 
form of various soliton expressions, we find that of eq. (|5.4[ ) may be written as the 
following determinant 



u 



= det 



1-t 2 (5-6) 
M m) L 4 = - ^? +m -, iJ = 1,2, ...,m 

L J 1 — XiXj 

as proved in Appendix A below. 

Eq. (|5.4|) provides the general solution to (|3.3|) a priori valid for any n (not necessarily 
positive or zero), that satisfies the convergence condition (if). Further imposing the initial 
conditions (i) also implies, by writing ( |3.3|) for n = —1, —2, ... — (m — 1) that i?_ m _i and 
i?_ m _ 2 diverge and that _R_ m _3, i?_ m _4, R-(2m-i) vanish, which in turn entail that 
= 0, = 0, •••' u _(2m-i) = ^' These latter conditions entirely fix the values of Aj 

to be 

Ai = II— — i = l,2,...,m (5.7) 



31^i J 
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At these values of the A's, the solution R n takes the following particularly simple form: 

w li---T W m)U n -{-3{Wi,...,W rn ) 

tin — -K77 7 TTT 7 r (,0-oj 

U n +\{Wl, Wm)U n + 2 (w 1 , ...,W m ) 

where 

U n (wi, w m ) = det [C7 n+2 j-2(M'i)]i<i,j< m ( 5 - 9 ) 
in terms of Chebyshev polynomials of the first kind expressed at values Wi = ^fxl + 1 / y/xi, 



% = 1,2, ...,to. Eqs.(|5.8| ) - (|5.9|) are proved in Appendix B below, by showing that oc 



U n (wi, io m ) where the proportionality factors drop out of the ratio giving R n . The 
initial conditions u^™) = for r = 1, 2, 2m — 1 translate for the choice (15.71) of the A's 



into the equivalent conditions that U- r (wi, ...,w m ) = 0. The latter conditions are clearly 
satisfied, as direct consequences of the reflection symmetry of the Chebyshev polynomials, 
namely that U- n -2(w) = —U n (w) for all n (in particular U-i(w) = 0), and therefore, as 
a determinant, U- r (wi, ...,w m ) always has at least two opposite lines (for even r) or a 
vanishing one (for odd r), provided r stays in the range 1, 2, .., 2m — 1. More precisely, this 
observation leads to a rank Int[(l + r)/2] for the matrix with determinant U- r (wi, ...,w r 
This gives an a posteriori justification of the choice (|5.7| ) for the A's. 



J m } 



6. Application: the hard dimer model 

6.1. The hard dimer model 

As a first application of the exact solution of Sect. 5, let us consider the so-called hard 
dimer model on tetravalent planar maps ||14j| . In this model, we enumerate tetravalent 
maps whose edges may be occupied or not by a dimer, with the hardness condition that 
no two adjacent edges may be simultaneously occupied by dimers. Assigning as usual a 
weight g per vertex and an activity z per dimer, the problem may be reformulated as that 
of enumerating maps with both tetra and hexa-valent vertices with a weight gi = 9 per 
tetravalent vertex and g% = 3zg 2 per hexavalent one. This is readily seen by shrinking 
the occupied edges, thus forming hexavalent vertices out of adjacent pairs of tetravalent 
ones, as shown in Fig. 6. Conversely, any given hexavalent vertex may be split into a pair 
of adjacent tetravalent ones in three different manners, hence the factor of 3. This is a 
particular case of the general model of previous section with m = 2. Note that in the hard 
dimer language, the notion of geodesic distance must be understood modulo the restriction 
that the minimal path between the two in- and out-coming legs must not cross any edge 
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(a) 



(b) 



-4 1 




4 1 





Fig.6: Example (a) of tetravalent two-leg diagram with hard dimers, repre- 
sented by thickened edges. The corresponding diagram obtained by shrinking 
the dimers (b) has both tetravalent and hexavalent vertices. The correspon- 
dence is three-to-one per dimer, as shown. 

occupied by a dimer. Indeed, after shrinking occupied edges it would correspond to going 
through an hexavalent vertex, which is forbidden in the language of graphs with tetra- and 
hexa-valent vertices. 

In the hard dimer case, eq.( [2.1|) reads 



2 D 3 



R=l + 3gR z + 30zg'R 



(6.1) 



while the recursion (13.31) reads 



R n — 1 + gR n (R n -i + R n + R n +i) + 3zg R n (R n - 2 R n -i + R n -i + Rn-iRn+i 

+ Rn+l + Rn+lRn+2 + Rn(2Rn-l + Rn + 2Rn+l)) 

The solution now reads 

U n (w ll w 2 )U n+ s(w ll w 2 ) 



R„ — R 



U n+1 (w 1 ,w 2 )U n+2 (w 1 , w 2 ) 



where 



U n {wi, w 2 ) = U n {wi)U n+2 {w 2 ) - U n+2 (wx)U n {w 2 ) 
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(6.2) 



(6.3) 



(6.4) 



and u>i = yfxl + 1/y/xi, with xi,x% the solutions with modulus less than one of the 
characteristic equation 



X2(x) = X2(x) = 1- gR[x+- +4) -3zg 2 R 2 (x 2 + 6x + 16 + - + (6.5) 



XI \ XX 2 



With this we find the initial terms 

„ w l + (x 1 + ^)(x 2 + j r ) l-AgR-A5zg 2 R 2 

ixn — -ft ; -1 ; ; r~ 7" — -rt" 



Rl — R 



l+(l + Xl + ±){l + x 2 + ±) ~ l-3gR-30zg 2 R 2 

1 - 7gR + (11 - 78z)g 2 R 2 + 23Azg 3 R 3 + \2\hz 2 g 4 R 4 
(1 - 3gR - 30zg 2 R 2 )(l - AgR - A5zg 2 R 2 ) 

1 - 8gR + (15 - 93z)g 2 R 2 + 366zg 3 R 3 + 1836z 2 g 4 R 4 



(1 - AgR - Ahzg 2 R 2 ){l - 7gR + (11 - 78z)g 2 R 2 + 23Azg 3 R 3 + 121hz 2 g 4 R 4 ) 

(6.6) 

Note that using eq. (|6.1| ) we may recast Rq = R—gR 3 (l + 15zgR) which is in agreement with 



the general formula ( |2.2|) for the two-leg diagram generating function, hence corroborating 
the relation Rq = r 2 - Finally, we may expand the first few GVs as 



G (g; z) = 1 + 2g + 3(3 + 5z)g 2 + 18(3 + 10z)g 3 + 18(21 + 105;? + 50z 2 )g 4 + ■■■ 
G^g; z)=g + A(2 + 3z)g 2 + (65 + 201z)g 3 + 2(277 + 1305^ + 585z 2 )g 4 + • • • (6.7) 
G 2 (g; z) = (l + 3z)g 2 + 3(5 + 21z)g 3 + (179 + 1005z + 531z 2 )g 4 + ■■■ 

As an illustration, we have represented in Fig. 7 the two-leg diagrams with two vertices and 
one hard dimer, corresponding to the terms of order g 2 z in eq.([ 



6.2. Tricritical point and fractal dimension 

The hard dimer model possesses a tricritical point characterized as follows. Introduc- 
ing the notation V = gR, the eq.( |6.1|) reads 



g = W(V) = V(l - 3V - 30zV 2 ) (6.8) 

The critical line of the model is characterized by the requirement that W'(V) = 0, leading 
to the parametric curve 

g = V ( 2 --v) 2 =i^ (6.9) 
y \3 J 90V 2 K 1 
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(b) 




(c) 

Fig. 7: The tetravalent two-leg diagrams with two inner vertices and one 
dimer are represented according to the distance n between the two legs, re- 
spectively n = (a), n = 1 (b) and n = 2 (c). Note that the two first 
diagrams of (c) have n = 2 as we forbid crossing the dimer. 

for < V < 1/3. The curve ends at a higher order tricritical point where, in addition we 
have W"(V) = 0, which fixes the critical values 



V, 



1 



1 



9c 



Rr 



(6.10) 



3 9 10 

We may approach this point by first setting z = z c = —1/10 and letting g approach 
g c = 1/9, in which case eq.( |6.8|) simplifies into 

3 N 



9 = 9c[l- 



Vc-V 

Vc, 



(6.11) 



We therefore set 



9 = 9c(l ~ e 



V = VJl-e 2 



(6.12) 
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and let e — > 0. At z = z c , the characteristic equation ( |6.5|) takes a particularly simple form 
when expressed in the variable y = x + 1/x — 2, namely 

X2(x) = = ^{Vyf - (1 - 3V)Vy + (1 - 3V) 2 (6.13) 

We see that at the tricritical point V = V c = 1/3, y = is a double root of eq.( |6.13p , hence 
x = 1 is a quadruple root of the characteristic equation. Moreover, eq.( |6.13|) amounts to 
the vanishing of a polynomial of degree 2 in the reduced variable u = yV/ (V c — V), namely 
P2(u) = u 2 /30 — u/3 + 1 = with the two roots u± = 5 ± i\fh. When approaching the 
tricritical point via ( |6.12j ), we therefore have the solutions 

y± = u±- o ( 6 - 14 ) 

1 — e z 

which imply the following leading behavior for the two x's with modulus less than 1 solving 
the characteristic equation: 



Xl2 = x ± = e~ ea± + 0(e 2 ) a± = \/5±zv / 5 (6.15) 

where the square roots are taken with positive real parts. As before, a sensible scaling 
limit is obtained by writing n = r/e in eq.( |6.3j ). This amounts to introducing a correlation 
length £ = e _1 = ((g c — g)/g c )~ 1 ^ 6 , which leads to a critical exponent v = 1/6. We 
may still interpret the quantity dp = 1/v = 6 as a tricritical fractal dimension as follows. 
Like in Sect. 4. 2, we consider the limiting ratio B n of (fOjj) where R n ,N now denotes the 
partition function for two-leg diagrams with vertices, distance at most n between the 
two legs and with hard particles with fugacity z = z c . A straightforward though tedious 
calculation yields 

1 (w + l)(n + 6) 
n 16632 (n + 3)(n + 4) 

x (33264 + 74088n + 63014n 2 + 26985?i 3 + 6215n 4 + 735n 5 + 35n 6 ) (6.16) 

5 6 
~ n 

2376 

in which we read off dp = 6. 



Further, expanding eq. (|6.3|) at order 2 in e, we obtain the tricritical scaling function 
T{r) counting surfaces with two marked points at a geodesic distance larger or equal to r: 

R — Rr, d 2 



Hr) = Ssa = ~ 2 d^ LogW ( sinh ( a+ D ' sinh ( a -£ 

(a 2 + - a 2 _) (a\ sinh 2 (a_|) - a 2 _ sinh 2 (a + §)) ^ 17 ) 
2 (a_ cosh (a~§) sinh — a+ cosh (a+§) sinh (a_|)) 2 
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where W(/i, fa) stands for the Wronskian determinant fify — faf[, and a± as in eq.( |6.15|) . 
As before, the scaling function Q(r) counting surfaces with two marked points at geodesic 
distance r is obtained by differentiating — J-(r). 

So far we have considered a particular approach to the tricritical point, namely with a 
fixed z = z c . Let us now show that the same exponent v is obtained for any other generic 
approach to the tricritical point and that moreover the same scaling function is obtained, 
up to a multiplicative redefinition of r. Indeed, expanding the relation g = W(V) around 
the critical values as g = g c + 5g, V = V c + 5V, z = z c + Sz, we obtain at leading order: 

8g + -£-6z = 3(5V) S (6.18) 
9 

Taking a generic approach 5g = —g c ae 6 , 5z = —z c (3e 6 , we arrive at 5V = — V c (a — (3)se 2 . 



Accordingly, the characteristic equation (|6.5| ) reads 

X2 ( X ) = = 3(5V) 2 P (-^) + 0(e d ) (6.19) 

which yields the same solution y± as in eq. ( |6.14j ) up to a multiplicative factor of (a — (3) s , 
which in turn leads to the same scaling function ( |6.17| ) up to a redefinition r — > r/{a — (3)^. 

One notable exception is when a = /3, or more precisely when z c Sg = g c Sz in which 
case the lhs of Q6.18 ) vanishes at leading order, and we must also consider higher terms 
in the expansion, namely of the form 5z5V. Such a situation occurs in particular if 



we approach the tricritical point along the critical line (|6.9|) . In that case, we have 
5g = — (SV) 2 . Now the appropriate scaling is 8g = — g c e 4 , and SV = —V c e 2 . The ex- 
pansion of the characteristic equation now leads to another polynomial equation of the 
form Qz(— yV c /SV)) = 0, where Q2(u) = u 2 /30 — u/3. This gives the solutions x\ = 1 and 
X2 = e _ae , with now a = y/lO. The appropriate scaling of the distance therefore still reads 
n = r/e, and we have an exponent v = 1/4 from the behavior of the correlation length 
£ = e _1 = ((<7c — g)/g c )~ 1 ^- This particular approach leads to a different scaling function 
T c which can be obtained from the form ( 6.17] ) by letting a+ — > and a_ = a = vTO, 
namely 

d 2 

F c (r) = "2^Log (v/fr cosh (^f r) - sinh (^f r)) (6.20) 

To conclude this section, we note that the generic tricritical scaling function T{r) of 
eq. ( |6.17| ) may be alternatively obtained by deriving the continuum counterpart of eq. ( |6.2p 
yielding the following differential equation for F(r) 

T {4 \r)-10J r (r)J r ''(r)-10J r ''(r)-5(J r '(r)) 2 + 10(J r (r)) 3 + 30(J r (r)) 2 + 30J r (r) = (6.21) 
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The scaling function J-(r) of eq. (|6.17|) is the unique solution to eq.( |6.21|) such that 
linir^o ^ 2 ^ r (r) = 6 while lim^oo J 7 if) = 0. The condition at r —>■ may be seen by 
writing T{r) = —2d 2 /dr 2 LogU(r), where U{r) is the continuum limit of u n and imposing 
the continuum counterpart of the discrete conditions u-i = U-2 = W-3 = 0, namely that 
U(0) = U'(0) = U"(0) = 0, hence U(r) oc r 3 , and F(r) ~ 6/r 2 . In the same spirit, the 
scaling function T c of eq. (|6.20| ) corresponding to a non-generic approach of the tricritical 



point along the critical line is the unique solution with the same boundary conditions as 
just mentioned of the truncated differential equation 

Tf\r) - 10^ c (r)^'(r) - 10^'(r) - 5(^(r)) 2 + 10(^ c (r)) 3 + 30(^ c (r)) 2 = (6.22) 

6.3. Continuum probability distribution for geodesic distances 

As in Sect. 4. 4, we may translate the above result into an explicit expression for the 
probability distribution for geodesic distances in tricritical surfaces with hard dimers. 
Again, we must compute the ratio R n ,N / Roo,n , where R n ,N now denotes the generat- 
ing function for two-leg diagrams with hard dimers at z = z c = —1/10, with inner 
vertices and legs distant by at most n. Starting from the contour integral ( f4.16| ) applied 
to the hard dimer problem, we perform the change of variables g = W(v), with W as in 
( |6.8|) , with the result 

/ dvW'(v) „ , , , . 

At large A, this integral is dominated by the saddle-point v = V c = 1/3. We now perform 
the change of variables v = V c (l — e _l7rsgn ^^/ 3 |u|/A 1 / 3 ) where u is to be integrated on the 
real line when A — > oo. This prescription is obtained by deforming the initial contour of 
integration around the origin of the v -plane so as to encompass only the image of the pole 
at v = 0. A suitable limit is obtained by further setting n = aA 1 / 6 , and using 



Rn(v) =R(l- e-^^T^ae-^Ju) + 0(4r)l 

V As iV~3 / 

= - fl - e- 1 ^ (1 + ^(ae-** + 0{-L)) 
9c V As As / 



(6.24) 



for u > and T as in eq.( |6.17|) and the corresponding counterpart for u < 0. We finally 

get 

R n ,N t duu 3 e- u Re(e-^(l + ^(ae^ v / ^))) (6.25) 

7rA3 Jo 
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This leads to the probability for tricritical two-leg diagrams to have a rescaled geodesic 
distance less than a: 

AT poo 

$( a )=6-y T -/ duu 3 e- u3 Re (e _i ^(l + ^(ae^v^))) (6-26) 

and by differentiating with respect to a, to the probability distribution p(pt) for distances 
equal to a. Note finally that for small a, we have $(a) ~ 5a 6 /1188, in agreement 
with eq. (|6.16| ), while the large a behavior obeys Fisher's law with a decay of the form 
exp(-Ca 6 / 5 ). 



7. General multicritical points: fractal dimensions and scaling functions 

7.1. Fractal dimensions 

In this section we consider the case of arbitrary m, in which case we may reach an 
m-th order multicritical point. For simplicity, we set g± = and we introduce the following 
notations 

92 = 9 9k = 9 k ~ X Zk (with z 2 = 1) V = gR 

k=2 V 7 

in terms of which eq.( |2.1|) reads simply W(V) = 9. The m-th order multicritical point is 
obtained by setting W = W" = ... = = which implies 



1 (m+l\ / a \ k-1 

2 W = (-D^^Li ( 1) (7. 2) 



III 



With these, eq. (|2.1|) simply reads 



g = w(v)= V^-^-p m+1 V c= ™ (7.3) 

while g c = W(V C ) = m/(6(m + 1)). For simplicity, we decide to approach the m-th order 
multicritical point by first setting all the z^'s to their critical values of eq.(|7.2|), and 
then letting g approach g c . As already shown in the previous section, this approach is 
sufficient to capture the generic multicritical behavior. As readily seen from eq.( [7.3| ), we 
must set 

g = g c (l - e 2m + 2 ) V = V c (l - e 2 ) (7.4) 
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The characteristic equation (5^2) takes a particularly simple form when expressed in 



(7.5) 



the variable y = x + 1/x — 2, namely 

This is a consequence of the following identity for Chebyshev polynomials 

A/2fc + l\ TT , . (2fc + l)!^ Z! ?/ .„ . 

where as before w = y/x + l/y/x. This identity is proved by showing that both sides satisfy 
the same recursion relation (3k+i — a 2 (3k = ( 2 fc fc j f 1 2 ), and have the same initial term /3q = 1. 
Eq.( [7T5"l) follows by substituting the identity (|7.6|) into eq.(|5.2|). With the value ( |7.3|) of 
W(V), the characteristic equation becomes 



Xm(x) = 0=1 — P, 



Pm(u) = 



Vr. J " m V Vr. ~ V 



,_ ~ U) (2/ + l)!(m-/)! 



(7.7) 



Approaching the multicritical point as in eq. (|7.4|) , we get the solutions 

e 2 

Vi=Ui- g (7.8) 

1 — e A 

where the Ui, i = 1, 2, ...,m are the distinct roots of P m . In turn, this yields the iCj's with 
modulus less than one 

Xi = e~ eai + 0(e 2 ) a, = y/ul (7.9) 

where the square roots are taken with positive real parts. Again, to get a sensible limit, 
we must write n = r/e in eq.( f5.8| ), which displays again a correlation length £ = e _1 = 
((9c — 9)l9c]~ v 'j with = l/(2m + 2). This leads to a generalized multicritical fractal 
dimension dp = 2(m + 2). 

7.1!. Scaling functions 

The expansion of eq. ( |5.8|) at order 2 in e gives the general multicritical scaling function 

jF(r) = -2-^LogW (sink (<H~) , sinh (a 2 , ...,sinh («m0) (7.10) 
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where W(/i, fa, f m ) is the Wronskian determinant det [f^ 1 ' > ] 1<ij<m - This scaling 
function may alternatively be obtained as the unique solution such that lim r ^o r 2 T{r) = 
2(2m — 1) and lim^oo T{r) = of the differential equation 

K m+1 [l + F] =U m+1 [l] (7.11) 



where 72. p [u] are the residues of the KdV hierarchy []T5[| , defined recursively by 7£' +1 = 
K™/4 - WR' p - u'Kp/2, n = 1/2. 

To show eq. (|7.11|) , let us recall that our recursion relation (|3.3|) coincides with that 
of the one-matrix model with even potential, up to a change 1 <-> n/Af. More precisely, 
writing formally eq. (|3.3| ) as ^ n ({Rj}; {gk}) = 1, the corresponding equation in the matrix 
model reads ^> n ({Rj}; {gk}) = n/Af. Approaching the m-th order critical point with 
the <7fc's given by eqs.( [7.1|) and (|7.2|) , and making the scaling ansatz g = g c (l — e 2m + 2 ) 5 
R n = -R c (l — e 2 u(ne)), the functional \E' n may be expanded at small e as 

VndRjh {9k}) = 1 + e 2m+2 (l - <Wu]) + 0(e 2m + 4 ) (7.12) 

where «S m +i[u] is a functional of u to be identified. In the context of the matrix model's 
double-scaling limit, we write ^ n ({Rj}; {gk}) = n/Af = 1 + e 2m + 2 (l — t), where t is the 
renormalized cosmological constant, which leads to the equation «S m +i[u] = t, where u is 
interpreted as the renormalized string susceptibility. From this, we identify <S m +i oc TZ m +i, 
as is well known from random matrix theory 0. Going back to our problem, the scaling 
ansatz Q4.25Q with R = R c (l — e 2 ) corresponds to having u = 1 + T and the equation 
ty n ({Rj}; {gk}) = 1 now yields <S m +i[l + T\ = 1. Eq.( [7.11|) follows by proportionality as 
we know from eq.( |2.1|) that T = must be solution. 

As an independent check, we may derive a characteristic equation 7r m (o 2 ) = by lin- 
earizing the continuum equation ( [7.11|) and looking for solutions of the form T{r) = e~ ar . 
Using the explicit value TZ P [1] = (-1)p( 2p )/2 2 p +1 , and the recursion relation for the KdV 
residues, we find the recursion relation 7r p (a 2 ) = (a 2 — 4)7r p _i(a 2 )/4 — (— l) p ( 2 p p )/2 2p+2 , 
with the initial value tyq = —1/4. Solving this recursion explicitly, we simply find that 
7r m ( a 2 ) = (2m + l)(-l) m " 1 ( 2 ^)/2 2m+2 P m (a 2 ), with P m as in eq.Q, hence the charac- 
teristic equation coincides with Xm( x ) = 0forx + l/x — 2 = e 2 a 2 . 

Two final remarks are in order. First, the scaling function (|7.10| ) is the generic scal- 



ing function for the m-th order multicritical point. It may also be obtained by starting 
with a problem involving vertices up to valences 2m' + 2, ml > m, and approaching any 
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point of the m-th order critical manifold of dimension m' — m, for which the characteristic 
equation has x = 1 as a 2m-times degenerate zero. Conversely, we may approach the 
m-th order multicritical point by first going to an m'-th order multicritical manifold with 
m' < m (which amounts to having x = 1 as a 2m/-times degenerate zero of the char- 
acteristic equation), and then approaching the m-th order multicritical point along this 
manifold. The corresponding scaling function is different from the generic one for the m-th 
order multicritical point and may be obtained as a limit of the form ( |7.10|) in which say 
ai, a m ' — > while aw+i, a m tend to the specific values coming from the restriction 
of the characteristic equation to the m'-critical manifold. 

7.3. Probability distribution for geodesic distances 

As in Sects. 4. 4 and 6.3, we may derive expressions for the probability distribution 
of geodesic distances in multicritical surfaces. Again, we compute the ratio R U) n I R<x,N, 

(c) 

where R n ,N denotes the generating function for two-leg diagrams at Zk = z k as given by 



eq.( [7.2| ), with iV inner vertices and legs at distance at most n. In the the contour integral 
( |4.16|) for the present multicritical problem, we perform the change of variables g = W(v), 
with W as in eq.( [7.1|) , to get, as before 



R 



n,N 



f dvW'iv) . , . . 



At large N, the integral is dominated by the saddle-point v = V c = m/6 and we perform 

■ 7rsgii(n) 1 

the change of variables v = V c (l — e % m + 1 \u\/N m + 1 ) where u is to be integrated on the 



i 



real line when iV — > oo. We obtain a suitable limit by further setting n = aN 2 ( m +^ , and 
using 



Rn(v) =R[l-e~ i ^ U 1 T(ae i ^+^ y /u)+Q( 1 2 )) 

g c V A^^+i v ' iV^+i 



(7.14) 



for u > with T as in eq.( [7.10| ) and the corresponding counterpart for u < 0. We finally 
get 

R n , N ~ {6m ^ )N ^ + 1)2 duu^e-- m+1 Re(e-^(l+H^e^^))) 

Tj-jVm + l JO ^ ' 

(7.15) 
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This gives the probability distribution for rescaled geodesic distances at most a 
$( a ) = (^+1) 2 [°° duu m+1 e- um+1 Re (V*WT (1 + ^(ae l rai v^)) 

cos (^i))r(^ T ) io v 

V 2(m+l) J \ m+1 ' 

(7.16) 

and that for rescaled geodesic distances equal to a by differentiation. As expected, we 
have <fr(a) oc a 2 ^ m+1 ^ at small a while p(a) decays for large a as exp(— Ca 2 ( m+1 )/( 27n+1 )). 



8. Discussion and conclusion 

In this paper, we have found explicit expressions for partition functions of planar maps 
with arbitrary even valences and with two marked points at a specified geodesic distance. 
From these we have extracted various scaling functions and critical exponents correspond- 
ing to various approaches to the different critical and multicritical points. We also obtained 
various probability distributions for distances in maps of fixed but large number of vertices. 
These results rely on a purely combinatorial technique based on bijections with trees. 

8.1. Relation to matrix models and integrability 

The results obtained here cannot apparently be reached by the more standard matrix 
integral approach, which naturally deals with the genus of the graphs rather than with 
geodesic distances. However, we observe strong analogies between the two problems, as 
apparent from the structure of our main recursion relation. Even more strikingly, in the 
continuum limit of the m-critical model, we found a differential equation of the form 
7£ m +i[u] = 7£ m +i[l] with u = 1 + JF, to be compared with the generalized KdV equation 
7£ m _i_i[u] = const, t for the all-genus string susceptibility u(£) where t plays the role of 
a renormalized cosmological constant in the double-scaling limit 0. This also amounts 
to having replaced the usual string equation [P,Q] = 1, Q = d 2 — u, d = d/dt, and 
P a differential operator of degree 2m + 1, by simply [P, Q] = 0. This modification is 
what makes the problem exactly solvable both in the discrete and continuum versions. In 
particular, we note the appearance of soliton-like expressions for our solutions, which are 
the sign of an underlying integrable structure. 

We expect that this scheme also applies to the general case of a c(p, q) CFT coupled 
to 2D quantum gravity, namely that the fundamental two-point scaling functions of the 
geodesic distance are still determined by differential equations of the form [P, Q] = 0, 
where P and Q are differential operators of respective degrees p and q. 
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8.2. Other two-point functions 

The continuum m-multicritical two-point function that we have derived was obtained 
as the scaling limit of the exact discrete solution. We expect on general grounds that 
it corresponds to the two-point correlation of the most relevant operator (pi in the corre- 
sponding c(2, 2m+l) conformal field theory coupled to 2D quantum gravity. It is clear that 
many more quantities may be derived from the general solution (|5.4|) . In the language of 
quantum gravity, this includes in particular correlations of the form ((pi(0)<fij(r)) involving 
less relevant operators <pj, j = 2,3, ...,m as functions of the geodesic distance r. In the 
language of this paper, natural discrete candidates for such correlators are the generating 

(k) 

functions G n for two-leg diagrams with legs distant by n and such that the outcoming leg 

^ , (fa's 

is attached to a 2/c-valent vertex. In this paper we have considered G n = + J^Gn , 

(k) 

a particular linear combination of the G„ . In the continuum limit, this combination has 
a non-vanishing overlap with the most relevant operator (pi of the theory, hence we may 
interpret our result as computing (<pi(0)(pi(r)) oc Q{r). We expect that by considering 

( k) 

other fine-tuned linear combinations of the Gn , we will obtain higher order terms in e, a 
term of order e 2j+1 corresponding to a two-point function (<pi(0)(pj (r)). Indeed, the above 
linear combinations amount to having inserted some operator <pj at the vertex adjacent to 
the outcoming leg, while the other leg remains untouched, and therefore still carries <pi. To 
compute these two-point functions, we observe that = = 9k(Q 2k ~ 1 )n-i,n, 

with notations as in eqs.( |3.4| )-( ]3.5|) . These are the same polynomials of the Rj's as those 
appearing in the recursion relation for the one-matrix model, and we may therefore apply 
the known results on the double-scaling limit in relation with KdV flows to identify 

(0i(O)^(r)) {X ~72 i [l + ^ (8.1) 



for j — 1, 2, m, with JF solving eq.(7.11 



These correlators also allow for an explicit exploration of the general problem of fractal 
dimensions, a subtle question already addressed in Refs.p, fT6| . 

8. 3. Generalizations 

We may adapt the techniques of this paper so as to address similar questions in more 
involved graph enumeration problems. Indeed, the case of planar maps with vertices of 
even valence that we considered here is a particular case of a more general class of planar 
maps called p-constellations. These are bipartite planar maps with say black and white 
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(a) (b) 

Fig.8: Example (a) of a two-leg 3-constellation with geodesic distance 1, and 



(b) the corresponding blossom tree (see Ref. [0 for precise definitions and 
characterizations) . 

vertices, such that all black vertices are p-valent, while the white ones have valences kp, 
k = 1,2, 3.... For p = 2, the black (two-valent) vertices may be wiped out while the white 
ones have arbitrary even valences, which is the case studied in this paper. For general 
p > 3, we may still consider two-leg diagrams with a univalent black vertex (in-coming 
leg) and a univalent white one (out-coming leg), not necessarily adjacent to the same 
face (see Fig.8 for an illustration). The geodesic distance between the legs is the minimal 
length of paths connecting them, with the constraint that one may only cross edges with 
the white vertex on the right and therefore the black vertex on the left, when going from 
the face adjacent to the in-coming leg to that adjacent to the out-coming one. These 
two-leg diagrams are in one-to-one correspondence with suitable bipartite blossom trees as 
shown in Ref. JT!^] . Using this equivalence, and limiting ourselves to white vertex valences 
up to (m+l)p for some fixed m we have shown that the generating function R n for two- leg 
diagrams with legs distant by at most n satisfies a polynomial recursion relation, which 
allows to determine them exactly. Remarkably, we find that the generic form for R n is 
again very similar to that of ( |5.4j) , except for a simple shift in the indices of , namely 

u {N) u {N) 

Kn ~ K (N) (AT) (*- 2 ) 
U n+l U n+p 
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where R = lim^^oo R n is a solution of a polynomial equation of degree (m + l)(p — 1), 
is a iV-soliton solution to KP of the form ( |5.5|) with N = (m+l){p— 1) — 1, e 77 * = — £™ +7V Ai, 
ccj being the solutions with modulus less than 1 of the characteristic equation associated 
to the recursion relation (and of degree N in x + l/x), and finally pi = Xi +xf + ... +xf , 
while 5j = 1/xi + 1/x^ + ... + \jx\~ 1 . Again, the integration constants Aj are determined 
by the initial conditions U-i = 0, i = 1,2, (2iV — 1), which in turn fully specifies the 
solution R n . Apart from these complications, the strategy of the present paper may be 
applied to these more involved cases of p-constellations. In particular, we have checked 



that we recover the expected generating function 1?2 = -Ro previously obtained in Refs.[17 



19] . However, we do not expect any new critical behavior or scaling functions other than 
that already obtained in this paper for p = 2. 

More promising, the same techniques can be applied to models of graphs with matter 
degrees of freedom such as the Ising model, whose study is under way ||20|| . and where new 
scaling functions appear. 

8.4- Connection to ISE 

Let us finally discuss the relation between the present work and the other combina- 
torial approach of Ref.0 for the tetravalent case. In this reference, a different bijection 
between tetravalent rooted maps and labeled trees is used, where the labels precisely en- 
code the geodesic distance from the root. Representing both the underlying tree and its 
labels by two discrete ID random walks, the problem translates into a discrete version of 
the so-called Brownian snake, which makes the connection between geodesic distances on 
tetravalent planar maps and random variables of the Integrated SuperBrownian Excursion 
(ISE). In particular, the continuum distribution $(«) for the simple critical case m = 1 
coincides with the distribution for the maximum of the ISE, whose suitably defined Laplace 



transform was obtained in Ref. |21| and coincides with T . Beyond this m = 1 case, it would 
be interesting to understand in ISE terms the counterpart of our multicritical results for 
arbitrary m. In particular, eq. (|7.16| ) may be inverted as 



T(a) = I" 4^e"» (l - * ( ) (8.3) 

(m + l)r(^ T )^o u~ V 

which displays the suitably generalized Laplace transform to be used in this case. In the 
ISE language, this suggests to consider excursions with "duration" u distributed according 

m + 2 

to a suitably generalized Ito law P{u) oc l/-u m + 1 . An example of the discrete version of 
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such "multicritical excursions" was considered in the context of (1 + 1) -D Lorentzian gravity 



in Ref. p2[ in the form of walks with steps of heights 1,2, • - • ,m with respective weights 
zi, ■ ■ ■ , z m which must be fine-tuned to some multicritical values. Understanding such a 
connection would provide a direct link between the so-called Lorentzian and Euclidean 
two-dimensional quantum gravities. 
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Appendix A. Proof of the determinant form ( |5.6|) for Un of eq.( |5.4| ). 

Eq. (|5.6|) is proved by expanding the determinant as 

1-x? 



ill / 

det Aff*> = s ^ II (<W) - X ^ +m j 



= £ S gn W £ y, n *i*s> n (-^< +m r^?-:> 

1=0 l<ii<...<ii<mj<£{i 1 ,...,i l } ie{ii,...,ii} 



1=0 l<i 1 <...<i i <m \r=l 







) det 




1 %i s %i t _ 





Ks,t<l 



(A.l) 

The identification with the form ( |5.4[ ) for is completed by computing the determinant 
in the last line of (|A.1| ) as 



det 



i-y 2 s 

1 - VsVt 



= n i 

l<s,t<l l< s <t<Z v 



y s ~ Vt 



VsVt 



(A.2) 



where y s = Xi a , obtained by direct application of the Cauchy determinant formula 

I [ (y s - yt)(z s - z t ) 



det 



upon taking z s = —l/y s . 



1 



y s + z t 



Ks<t<l 



Ks,t<l 



(A.3) 



Ks,t<l 
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Appendix B. Proof of the determinant form ( |5.8|) -( |5T9|) . 

Let us start from the determinant formula ( |5lf ) in which we substitute the values 
(O) of the A;: 

M M = det M^ m) 

|3 1 - XiXi 

Mfl, < = Sij - *™ +m ^ /'•./ = 1.2 m 



(B.l) 



1 h3 



J]_Xi- Xl 

We now multiply M^™' to the right by the matrix _p( m ) with entries 

r -i ^ m -j 

pirn) 



X- 



-■'■■) Ylxi-xi 



(B.2) 



This gives the general term 



X 



m-j 



k=l l^k 



1 - XiXi 



2LL 



(B.3) 



Defining 



1 — XX\ 



(B.4) 

k=l l^k x k 

we see that <p(x) is the Lagrange interpolation polynomial of degree less or equal to m—1, 
that takes the values cpj(l/xk) = x^~ 3 at the m distinct points x = l/x^, for k = 1, 2, m. 
We conclude that <Pj(x) = x 3 ~ for j = 1, 2, m. This allows to finally rewrite ( |B.3| ) as 



M {m) p {m) 



X- 



m+ - 



« n 



I o 

a* - ^z Vx, 2 



n — 1 , ■ 

— +o 



m+ - 



(B.5) 



X: 



\fx~i U n+2 j-2(Wi) 



Y]_Xi-Xi 

with Wi = ^/xl+ 1/y/xi as before. Taking the determinant, we arrive at 



u 



(m) 



1 



det P( m ) 



n 



i=l 



\ 



V 



x det [I7 ri+ 2j-2(w»)] 1 < i 



(B.6) 



/ 



which immediately leads to eqs.(|5.8|)-(p\9|), by noting that all the prefactors drop out of 
the ratio. 
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